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Abstract. ROSAT HRI observations have been used to determine an upper limit of the Crab 
pulsar surface temperature from the off-pulse count rate. For a neutron star mass of 1.4 and 
a radius of 10 km as well as the standard distance and interstellar column density, the redshifted 
temperature upper limit is T s °° < 1.55 X 10 6 K (3cr). This is the lowest temperature upper limit 
obtained for the Crab pulsar so far. Slightly different values for T s °° are computed for the various 
neutron star models available in the literature, reflecting the difference in the equation of state. 

1. Introduction 

Associated with the historical supernova of AD 1054, the Crab pulsar PSR 0531+21 
is the youngest neutron star of a population of over 550 known rotation powered 
pulsars, and the only one, for which the age is known with high accuracy. The ques- 
tion about the thermal properties of the Crab pulsar has therefore been of special 
importance for the pulsar and neutron star theories since its discovery; knowing the 
surface temperature or temperature upper limit allows a comparison with the pre- 
dictions of cooling theories for a neutron star at an early stage of its evolution and 
may yield important constraints for the equations of state and the internal struc- 
ture of the neutron star. However, although the Crab is one of the most frequently 
observed sources in high energy astrophysics, previous experiments have failed to 
detect any spectral evidence for thermal emission from the neutron star's surface. 
The intense magnetospheric emission during the pulse-on phase prevents the detec- 
tion of the thermal radiation associated with the neutron star's initial heat content, 
and during the pulse-off phase the thermal emission from the neutron star surface 
is buried under the bright synchrotron emission of the surrounding nebula. Only 
upper limits of the Crab pulsar surface temperature T s could be determined for that 
reason. First results in this respect were reported by Wolff et al (1975) and Toor 
& Seward (1977). Based on rocket born lunar occultation observations of the Crab 
supernova remnant in 1974, T s < 4.7 x 10 6 K and T s < 3 x 10 6 K, respectively, 
were inferred from the off-pulse count rate upper limit, which, interpreted in terms 
of blackbody radiation from the entire neutron star, is equivalent to a limit of the 
neutron star surface temperature. An improved result became available from the 
Einstein Observatory using the High Resolution Imager (HRI) few years later. As- 
suming a column density of 3 x 10 21 cm -2 and a distance of 2 kpc, a 3c upper limit 
of T s < 2.5 x 10 6 K was reported by Harnden & Seward (1984) for a neutron star 
of radius R=10 km and mass M=1.4 M @ . The lowest temperature upper limit re- 
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ported for the Crab pulsar so far, was derived from an analysis of the 1969 and 1975 
glitches and post-glitch behaviour and the interpretation in terms of unpinning and 
repinning of crustal superfluid vortex lines. The creep relaxation time obtained from 
a fit to the post-glitch timing data implies an internal temperature of ~ 3 x 10 8 K, 
which for a 1.4 M @ neutron star with a reasonable stiff equation of state results in 
a redshifted surface temperature upper limit of T s °° < 1.6 x 10 6 K (Alpar et al 1985). 

In the present work, ROSAT HRI data from the Crab SNR are used to revisit the 
surface temperature of PSR 0531+21. Since the performance of the ROSAT X-ray 
telescope and HRI exceeds that of the Einstein Observatory in terms of the HEW 
angular resolution and scattering (Aschenbach 1988, David et al 1992), the pulsar 
can be separated more clearly from the surrounding nebula emission. The result is 
an improved upper limit, which for the first time constraints the thermal evolution 
models of neutron stars. 

2. Observations and Data Analysis 

The Crab SNR was one of the first targets observed in the ROSAT pointing program 
between 20-24 March 1991. During 8,785 s effective exposure with the HRI in the fo- 
cus of the XRT, ~ 2.35 x 10 6 cts were collected within a 2 arcmin wide circle centred 
on the pulsar. About 95% of the total emission was found to come from the intense 
synchrotron nebula (cf. Rosat Calendar 12/1992). To determine the off-pulse count 
rate from the DC level of the Crab pulsar's light curve, 111,366 counts were selected 
from a 6 arcsec wide aperture around the pulsar. After barycentring and correcting 
the photon arrival times using the approach described in Becker et al (1993a), each 
X-ray photon arrival time was related to the pulsar phase <f>, using the Crab Pulsar 
Timing ephemeris of 15. March 1992 (Lyne & Pritchard 1992). The corresponding 
Crab pulsar light curve in the ROSAT energy band is shown below. 



Fig.l Folded X-ray light 
curve of the Crab pul- 
sar from ROSAT HRI ob- 
servations of 20-24 March 
1991. A single cycle 
is plotted in 333 phase 
bins of ~ 100 /xsec width. 
The DC level is fitted to 
(84.5 ± 1.2) cts/bin, indi- 
cated by the dashed line. 
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The (0.1-2.4) keV DC flux, integrated over the observation time, was found to be 
(84.5 ± 1.2) cts/bin by fitting a constant level to the off-pulse part of the X-ray light 
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curve, shown in figure 1. However, the ROSAT HRI point spread function contains 
only 72% of the total energy from a point source within 6 arcsec radius by which 
the DC flux has been corrected to (84.5 + 1.2) x 100/72 < 119 cts/bin. Since the 
X-ray image during the ofF-pulse phase does not show any evidence for a residual 
point-like source, the counts observed during this period are likely to be associated 
with the nebula rather than with the pulsar. Therefore, only an upper limit of 0.069 
cts/s (3a) for the Crab pulsar ofF-pulse flux can be derived (Becker 1994). 

3. Results 

Assuming blackbody emission from the surface of the neutron star of radius R, as 
well as the distance d and absorption column density Nh, the surface temperature 
T s can be derived from the observed count rate. However, due to the gravitational 
redshift z only T s °° = T s /(1 + z) can be observed. For R=10 km, M=\A M@, 
N H = 3 x 10 21 cm -2 and d=2 kpc the 3cr upper limit of 0.069 cts/s corresponds 
to a 3cr temperature upper limit of T s °° = 1.55 x 10 6 K or T s = 2 x 10 6 K. This 
result is in line with the upper limit derived by Alpar et al (1985) and is about 
0.5 x 10 6 K lower than the limit reported by Harnden & Seward (1984). The general 
relation between the redshifted surface temperature and the neutron star radius or 
bolometric luminosity L^ oV respectively, is shown in figure 2 for different interstellar 
column densities Nh- 




10 20 30 40 50 33.5 34.0 34.5 

Radius R (km) Log Lt 0l (erg/s) 



Fig. 2a Neutron star surface temperature T, 00 vs. i? = i?oo/(l + z), inferred 
from the Crab pulsar off-phase count rate upper limit for Nh = (li3,5) x 
10 21 atoms/cm . (b) Bolometric luminosity vs. = T,/(l + z). A distance 
of d = 2 kpc and M=1A M© have been assumed. 

In the literature various neutron star models of different equations of state have 
been proposed. To compare our result with the predictions from the corresponding 
cooling curves, the stellar parameters M and R have been used for computing T s 
and T s °°, respectively. Table 1 shows the results derived for the different models 
which are discussed, for example, in Nomoto & Tsuruta (1987), Van Riper (1991), 
Umeda et al 1993, Chong & Cheng (1993) and Page (1994). The errors quoted for 
T and L correspond to the uncertainty of d={ 2 ± 0.5 ) kpc and Nh = ( 3 ± 1 ) x 

91 / 2 

10 atoms/cm . 
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TABLE 1 3cr Temperature Upper Limits for PSR 0531+21 



Model Radius Mass T s x 10 6 T s °° x 10 6 log T s °° log 

km M@ K K K erg/s 



PS 1 


16 


10 


1 


31 


1 


61 


+0 
-0 


19 
21 


1 


40 


+0 
-0 


18 
16 


6 


15 


+0 
-0 


05 
06 


33 


97 


+0 
-0 


19 
24 


PS 2 


15 


83 


1 


40 


1 


63 


+0 
-0 


19 
21 


1 


40 


+0 
-0 


18 
16 


6 


15 


+0 
-0 


05 
06 


33 


97 


+0 
-0 


19 
24 


MPA 3 ' 4 


12 


45 


1 


40 


1 


81 


+0 
-0 


21 
24 


1 


48 


+0 
-0 


19 
17 


6 


17t°o 


05 
06 


33 


90 


+0 
-0 


19 
24 


PAL33 4 


11 


91 


1 


40 


1 


85 


+0 
-0 


22 
24 


1 


49 


+0 
-0 


19 
18 


6 


17 


+0 
-0 


05 
06 


33 


89 


+0 
-0 


19 
24 


UV14 4 ' 5 


11 


20 


1 


40 


1 


90 


+0 
-0 


23 
25 


1 


51 


+0 
-0 


20 
18 


6 


18 


+0 
-0 


05 
06 


33 


87 


+0 
-0 


19 
24 


UU 6 


11 


14 


1 


40 


1 


91 


+0 
-0 


23 
25 


1 


51 


+0 
-0 


20 
18 


6 


18 


+0 
-0 


05 
06 


33 


87 


+0 
-0 


19 
24 


PAL32 4 


11 


02 


1 


40 


1 


92 


+0 
-0 


23 
25 


1 


52 


+0 
-0 


20 
18 


6 


18 


+0 
-0 


05 
06 


33 


87t o 


19 
24 


FP 1 


10 


90 


1 


29 


1 


90 


+0 
-0 


23 
25 


1 


53 


+0 
-0 


20 
18 


6 


18 


+0 
-0 


05 
06 


33 


85 


+0 
-0 


20 
24 


pp4,5 


10 


85 


1 


40 


1 


93 


+0 
-0 


23 
25 


1 


52 


+0 
-0 


20 
18 


6 


18 


+0 
-0 


05 
06 


33 


86 


+0 
-0 


19 
24 


AV14 4 ' 5 


10 


60 


1 


40 


1 


96 


+0 
-0 


23 
26 


1 


53 


+0 
-0 


20 
18 


6 


18 


+0 
-0 


05 
06 


33 


86 


+0 
-0 


20 
24 


AU 6 


10 


40 


1 


40 


1 


98 


+0 
-0 


24 
26 


1 


53 


+0 
-0 


20 
18 


6 


19 


+0 
-0 


05 
06 


33 


85 


+0 
-0 


20 
24 


FP(pion) 1 


9 


40 


1 


27 


2 


04 


+0 
-0 


25 
27 


1 


58 


+0 
-0 


21 
19 


6 


20 


+0 
-0 


05 
06 


33 


82 


+0 
-0 


20 
24 


PAL(kaon) 7 


8 


20 


1 


40 


2 


27 


+0 
-0 


27 
30 


1 


60 


+0 
-0 


21 
19 


6 


20 


+0 
-0 


05 
06 


33 


80 


+0 
-0 


20 
24 


BPS 1 


7 


90 


1 


35 


2 


30 


+0 
-0 


28 
30 


1 


62 


+0 
-0 


21 

20 


6 


21 


+0 
-0 


05 
06 


33 


79 


+0 
-0 


20 
24 


BPS 2 


7 


35 


1 


40 


2 


45 


+0 
-0 


30 
32 


1 


62 


+0 
-0 


21 

20 


6 


21 


+0 
-0 


05 
06 


33 


78 


+0 
-0 


20 
25 


PAL(kaon) 7 


7 


00 


1 


40 


2 


55 


+0 
-0 


31 
34 


1 


63 


+0 
-0 


21 

20 


6 


21 


+0 
-0 


05 
06 


33 


78 


+0 
-0 


20 
25 



The stellar parameters Radius and gravitational Mass are taken from: 1 Umeda et al 
(1993), 2 Van Riper (1991), 3 Miither et al (1987), 4 Page (1994), 5 Wiringa et al (1988), 
6 Chong & Cheng (1993), 7 Thorsson et al (1994). 



The comparison between the predicted cooling curves of the PS model and the 
surface temperature upper limit derived for the Crab pulsar is shown in figure 3. 
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Fig. 3 3<j surface temperature up- 
per limit for PSR 0531+21, using 
the PS model parameters. Also 
shown are the cooling curves from 
Umeda et al (1993). Dotted and 
dashed lines represent the ther- 
mal evolution with internal fric- 
tional heating for strong (J3), weak 
(J2) and superweak (Jl) pinning of 
crustal superfiuid vortex lines. 
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4. Conclusion 

Due to the superior performance of the ROSAT telescope an improved upper limit 
of the Crab pulsar surface temperature has been obtained. Whereas the previous 
upper limit derived from data of the Einstein Observatory was consistent with the 
prediction of standard neutron star cooling, even if heating processes like frictional 
heating or crust cracking have been taken into account (Van Riper 1991, Page 1992, 
Umeda et al 1993, Chong & Cheng 1993). Even more, c.f. figure 3, the temperature 
upper limit obtained with ROSAT is well below the temperature predicted by the 
PS model which invalidates its application. A comparison with the FP model shows 
agreement with the observation but only if the uncertainty in the distance and in 
the column density is stretched to the limits. Strong frictional heating in the FP 
model is also excluded by our observation. 

A more detailed comparison between the different thermal evolution models and 
the Crab pulsar surface temperature upper limit is in preparation (Becker 1994). 
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